The quantum Ji-J[-J2 spin-1 Heisenberg model: Influence of the 
interchain coupling on the ground-state magnetic ordering in 2D 
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Abstract. - We study the phase diagram of the isotropic J1-J1-J2 Heisenberg model for spin-1 
particles on an anisotropic square lattice, using the coupled cluster method. We find no evidence 
for an intermediate phase between the Neel and stripe states, as compared with all previous 
results for the corresponding spin-1/2 case. However, we find a quantum tricritical point at 
J[/Ji ^ 0.66 ± 0.03, J2/J1 ~ 0.35 ± 0.02, where a line of second-order phase transitions between 
the quasi-classical Neel and stripe-ordered phases (for J[/Ji < 0.66) meets a line of first-order 
phase transitions between the same two states (for J[/Ji > 0.66). 



Over the last 20 years or so much theoretical effort [1-8] 
has been expended on the J1-J2 model in which the spins 
situated on the sites of a two-dimensional (2D) square 
lattice interact via competing isotropic Heisenberg in- 
teractions between the nearest-neighbour (Ji) and next- 
nearest-neighbour (J2) pairs. The exchange bonds Ji > 
promote antiferromagnetic order, while the J2 > bonds 
act to frustrate or compete with this order. Such frus- 
trated quantum magnets continue to be of interest because 
of the possible spin- liquid and other such novel phases that 
they can exhibit (see, e.g., Ref. [9]). 

The syntheses of compounds that can be closely de- 
scribed by the spin-1/2 version of the model, such as 
Li2VO(Si,Ge)04 [10] and VOM0O4 [11] have further fu- 
elled theoretical interest. It is now widely accepted 
that the spin-1/2 J1-J2 model on the 2D square lattice 
has a ground-state phase diagram showing two phases 
with quasi-classical long-range order (LRO) (viz., a Neel- 
ordered phase at small values of J2/J1 and a coUinear 
stripe-ordered phase at large values of J2/J1), separated 
by an intermediate quantum paramagnetic phase without 
magnetic LRO in the parameter regime < a < a^, 
where a = J2/J1 and a J « 0.4, ~ 0.6. Furthermore, it 
has been argued recently that the quantum phase transi- 
tion between the quasi-classical Neel phase and the quan- 
tum paramagmetic phase present in the 2D J1-J2 model 
is not described by a Ginzburg-Landau type critical the- 



ory, but rather may exhibit a deconfined quantum critical 
point [12]. Other authors [13] have argued that the transi- 
tion is not of this second-order type due to the deconfine- 
ment of the fractionally-charged spinous, but is rather a 
(weakly) first-order transition between the Neel phase and 
a valence-bond solid phase with columnar dimerisation. 

Such frustrated quantum magnets often have ground 
states that are macroscopically degenerate. This feature 
leads naturally to an increased sensitivity of the under- 
lying Hamiltonian to the presence of small perturbations. 
In particular, the presence of anisotropies in real systems 
that are well characterised by the J1-J2 model, either in 
spin space or in real space, naturally raises the issue of 
how robust are the properties of the J1-J2 model against 
any such perturbations. There have been several recent 
studies addressing this question. For example, in the case 
of spin anisotropies, generalizations of the J1-J2 model 
have been studied for the spin-1/2 case, in which either 
the frustrating next-nearest-neighbour interaction or the 
nearest-neighbour interaction is anisotropic [14,15]. 

In the alternative case of real-space anisotropies, for ex- 
ample, a recent study [16] investigated the effects of in- 
cluding an interlayer coupling (Jl) for the spin- half Ji- 
J2 model on a stacked square lattice. In a previous pa- 
per of our own [17] we moved instead in the direction of 
one-dimensionality by investigating a spin-half spatially 
inhomogeneous Ji-J[-J2 model in which the nearest- 
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neighbour bonds on the square lattice differ for the in- 

trachain (Ji) and interchain {J[) cases. The model can 
thus be viewed as parallel (Ji) chains, coupled by nearest- 
neighbour (J{) and next-nearest-neighbour (J2) interac- 
tions that frustrate each other. We found the surprising 
and novel result that for the spin- 1/2 case there exists 
a quantum triple point below which there is a second- 
order phase transition between the quasi-classical Nccl 
and stripe-ordered phases with magnetic LRO, whereas 
only above this point are these two phases separated by 
the intermediate magnetically disordered liquid-like phase 
seen in the pure spin-1/2 J1-J2 model (i.e., at J[ = J\). 
The quantum triple point was found to occur at J(/Ji ~ 
0.60 ± 0.03, J2/J1 « 0.33 ± 0.02. 

In the present work we extend the study of the J1-J1-J2 
model to consider the spin-1 case, which is computation- 
ally more challenging than the previous spin-1/2 case. As 
in the previous case we again use the much-studied coupled 
cluster method (CCM). Our main rationale for the present 
stiidy is that one knows in general that the spin quan- 
tum number can play an important and highly non-trivial 
role in these strongly correlated magnetic-lattice systems, 
which often exhibit rich and interesting phase scenarios 
due to the interplay between the quantum fluctuations 
and the competing interactions present. The strength of 
the quantum fluctuations can be tuned cither by intro- 
ducing spin-anisotropy terms in the Hamiltonian [18] or 
by varying the spin quantum number s [19]. 

While the general trend is that as s is increased the 
effects of quantum fluctuations reduce, one also knows 
that there can be significant deviations from it. A par- 
ticularly well-known example is the since-confirmed pre- 
diction of Haldane that integer-spin systems on the linear 
chain would have a nonzero excitation energy gap, whereas 
half-odd- integer spin systems would be gaplcss [20]. In- 
deed, such deviations from general trends provide one of 
the main reasons why quantum spin-lattice problems still 
maintain such an important role in the general investiga- 
tion of quantum phase transitions. 

In this context we note that the recent discovery of su- 
perconductivity with a transition temperature at Tc ~ 
26 K in the layered iron-based compound LaOFeAs, when 
doped by partial substitution of the oxygen atoms by flu- 
orine atoms [21], has been followed by the rapid discovery 
of superconductivity at even higher values of Tg (> 50 K) 
in a broad class of similar quaternary compounds. Enor- 
mous interest has thereby been engendered in this class 
of materials. Of particular relevance to the present work 
are the very recent first-principles calculations [22] show- 
ing that the undoped parent precursor material LaOFcAs 
is well described by the spin-1 J\-Ji model on the square 
lattice with J\ > 0, J2 > 0, and J2/J1 ~ 2. 

Returning to our present system, we note that while 
the s = 1/2 version of the J\-J'x-Ji model under discus- 
sion has been studied by various groups [17,23-26], very 
few calculations have been performed on the s = 1 case 
up till now. An exception is the two-step density-matrix 



renormalisation group study of Moukouri [26,27] that we 
discuss later in our concluding remarks. It has also been 
observed that quantum fluctuations can destabilize the or- 
dered classical ground state (GS), even for values s> 1/2, 
for large enough values of the frustration [1,28]. 

The model itself comprises a set of ^ 00 spin-1 par- 
ticles on a spatially anisotropic square lattice interacting 
via isotropic Heisenberg couplings, but with three kinds 
of exchange bonds. Its Hamiltonian is given by 

a = Ji Sj^i ■ Sj+i,; + J'l Sj^i ■ Sj^i+i 

id i,l 
i,l 

where the index (i, /) labels the x (row) and y (column) 

components of the lattice sites. The exchange constant 
Ji (which we henceforth set to 1) measures the intrachain 
bond strength along the row direction, while J( and J2 are 
the transverse (column) and diagonal interchain couplings 
respectively. The model retains the basic physics of the 
J1-J2 model (that is recovered when J[ = Ji), and has 
two types of classical GS, namely, the Neel (tt, tt) state 
and stripe states (columnar stripe (tt, 0) and row stripe 
(0, tt)). There is clearly a symmetry under the interchange 
of rows and columns, Ji ^ J[. which implies that we need 
only consider the range of parameters with J{ < Ji. The 
(first-order) classical phase transition occurs at the point 
of maximal frustration, J| = J[/2, VJi > J(. 

The CCM (see, e.g., Refs. [29-31] and references cited 
therein) employed here is one of the most powerful and 
most versatile modern techniques in quantum many-body 
theory. It has been successfully applied to various quan- 
tum magnets (sec Refs. [16, 31-35] and references cited 
therein). The CCM is particularly appropriate for study- 
ing frustrated systems, for which the main alternative 
methods are often only of limited usefulness. For example, 
quantum Monte Carlo techniques arc particularly plagued 
by the sign problem for such systems, and the exact diago- 
nalisation method is restricted in practice, particularly for 
s > 1/2, to such small lattices that it is often insensitive 
to the details of any subtle phase order present. 

We now briefly describe the CCM means to solve 
the ground-state (gs) Schrodinger ket and bra equations, 
i/]*) = i?]*) and = E{^>\ respectively (and see 

Refs. [29-34] for further details). The flrst step in im- 
plementing the CCM is always to choose a model state 
)$) on top of which to incorporate later in a systematic 
fashion the multispin correlations contained in the exact 
ground states j^") and {^\. More specifically, the CCM 
employs the exponential ansatz, j^E') =c"^|<I>) and (^P) = 
(<I']iSe~'^. The correlation operator S is expressed as 5 = 
^j_^QiS/C/ and its counterpart is 5* = 1 + '^Zi^a'^i^^T ■ 
The operators C/ = (Cj")^ with Cq = 1, have the prop- 
erty that ($|C+ = ; V/ ^ 0. They form a complete 
set of multispin creation operators with respect to the 
model state ]$). The ket- and bra-state correlation co- 
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efficients (5/, 5/) are calculated by requiring the gs en- 
ergy expectation value H = (^'|i?|\E') to be a minimum 
with respect to each of them. This immediately yields 
the coupled set of equations {^\CY He^ \^) = and 
i^lSie-'^He'^ - -E)C/|$) = 0; V/ 7^ 0, which we solve 
in practice for the correlation coefficients {Si,Si) within 
specific truncation schemes described below, by making 
use of parallel computing routines [36] . 

In order to treat each lattice site on an equal footing we 
perform a mathematical rotation of the local spin axes on 
each lattice site such that every spin of the model state 
aligns along its negative z-axis. Henceforth our description 
of the spins is given wholly in terms of these locally de- 
fined spin coordinate frames. In particular, the multispin 
creation operators may be written as C/ = .s^ • • • , 
in terms of the locally defined spin-raising operators = 
sf + on lattice sites i. Having solved for the multispin 
cluster correlation coefficients {Si,Si) as described above, 
we may then calculate the gs energy E from the relation 
E = {^\e~^ He^\^), and the gs staggered magnetisation 
M from the relation M = - X^^i sf |*) which holds 
in the rotated spin coordinates. 

Although the CCM formalism is clearly exact if a com- 
plete set of multispin configurations {/} with respect to 
the model state |$) is included in the calculation of the 
correlation operators S and S, in practice it is necessary 
to use systematic approximation schemes to truncate them 
to some finite subset. In our earlier paper on the s = 1/2 
version of the present model [17], we employed, as in our 
previous work [16,31-34], the locahsed LSUBn scheme in 
which all possible multi-spin-flip correlations over differ- 
ent locales on the lattice defined by n or fewer contiguous 
lattice sites are retained. 

However, we note that the number of fundamental 
LSUBn configurations for s ~ I becomes appreciably 
higher than for s = 1/2, since each spin on each site i 
can now be flipped twice by the spin-raising operator sf . 
Thus, for the s = 1 model it is more practical, but equally 
systematic, to use the alternative SUBn-m scheme, in 
which all correlations involving up to n spin flips span- 
ning a range of no more than m adjacent lattice sites are 
retained [31,37]. We then set m ~ n, and hence em- 
ploy the so-called SUBr? n schc^mc\ More; generally, the 
LSUBm scheme is thus equivalent to the SUBn m scheme 
for n ~ 2sm for particles of spin s. For s — 1/2, LSUBn = 
SUBn-n; whereas for s = 1, LSUBn = SUB2n-n. The 
numbers of such fundamental configurations (viz., those 
that are distinct under the symmetries of the Hamiltonian 
and of the model state [$)) that are retained for the Neel 
and stripe states of the current s = 1 model at various 
SUBn-n levels are shown in Table 1. 

Although we never need to perform any finite-size scal- 
ing, since all CCM approximations are automatically per- 
formed from the outset in the N ^ 00 limit, we do need 
as a last step to extrapolate to the n ^ 00 limit in the 
truncation index n. We use the same well-tested scaling 
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Table 1: Number of fundamental configurations (U f.c) for 
the SUBn-^ (n = {2, 4, 6, 8}) scheme for the Neel and stripe 
columnar states, for the spin-1 J\-J[-J2 model. 



Method 


Hf.c 


SUBn-n 


Neel 


stripe 


SUB2 2 


2 


1 


SUB4 4 


28 


21 


SUBG () 


744 


585 


SUB8-8 


35629 


29411 



-1.4 



-1.6 




-2.2 
-2.4 



Fig. 1: Extrapolated CCM SUBn-^r results for the gs energy 
per spin, E/N , of the spin-1 Ji-J{-J2 model (with Ji = 1), 
for J[ = 0.2, 0.4, 0.6, 0, 7, 0.8, 0.9, 1.0. The SUBn-n results are 
extrapolated in the limit n — > 00 using the set n = {2, 4, 6, 8}. 



laws as for the s = 1/2 model for the gs energy per spin 
E/N and the gs staggered magnetisation M, 

E/N = ao + ain"^ + a^n'^ , (2) 

M = 60 + (61 + 62n"') • (3) 

We report below on two separate sets of CCM calculations 
for this model, for given parameters ( Ji = 1, J(, J2), based 
respectively on the Neel state and the stripe state as the 
model state [<!>). In each case we have 4 calculated data 
points to fit the 3 unknown parameters in eqs. (2) and (3). 
Thus, we present below our final results ao iov E/N and 60 
for M, from fitting to the above schemes with the SUBn- 
n solutions for n = {2,4,6,8} as input. We note that 
for the corresponding s = 1/2 model we could perform 
LSUBn = SUBn-n approximation calculations for n = 
{2, 4, 6, 8, 10}. This enabled us to perform extrapolations 
using the sets n = {2, 4, 6, 8} and n = {2, 4, 6, 8, 10} as well 
as the preferred set n = {4, 6, 8, 10}. Gratifyingly, all sets 
yielded very similar extrapolated results, even near phase 
boundaries and the quantum triple point, which gave us 
great confidence in the accuracy and robustness of our 
extrapolation scheme. 
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Figure 1 shows the gs energy per spin as a function of J2 
for various vahics of J(, extrapolated from the raw CCM 
data as discussed above. The raw SUBn-n data termi- 
nate at some particular values. This occurs for the CCM 
curves based on both the Nccl state and the stripe state as 
the model state |$). Such a termination point arises due 
to the solutions of the CCM equations becoming complex 
at this point, beyond which there exist two branches of 
complex-conjugate solutions [31]. In the region where the 
solution reflecting the true physical situation is real, there 
actually also exists another real solution. However, only 
the (shown) upper branch of these two solutions reflects 
the true physical situation, whereas the lower branch does 
not. The physical branch is easily identified in practice as 
the one which becomes exact in some known (e.g., pertur- 
bative) limit. This physical branch then meets the cor- 
responding unphysical branch at some termination point 
beyond which no real solutions exist. The termination 
points shown in fig. 1 are the extrapolated n — > 00 ter- 
mination points and are evaluated using data only up to 
the highest level of the CCM approximation schemes used 
here, namely SUB8-8 for the s = 1 model. The SUBn-n 
termination points are also refiections of phase transitions 
in the real system, as we discuss more fully below. 

The maxima of the extrapolated gs energy curves are 
close to the corresponding classical transition points at 
J| = 0.5 J{. The extrapolated gs energy curves of the Neel 
and stripe states meet smoothly with the same slope for 
J'l < 0.66 ±0.03. This behaviour is indicative of a second- 
order phase transition. By contrast, for J[ > 0.66 ± 0.03 
the behaviour is typical of a first-order phase transition 
where the curves cross with a discontinuity in the slope. 
A comparison of fig. 1 for the present s = 1 model with 
fig. 2 of Ref. [17] for the s = 1/2 model clearly shows the 
distinct differences between the two cases. Thus, for the 
s = 1/2 case each gs energy curve for the Neel state for 
values J{ > 0.60 terminates before it can reach the cor- 
responding gs energy curve for the stripe state within the 
region that reflects the true physical situation (viz., where 
the calculated staggered magnetisation is positive), indi- 
cating the opening up of an intermediate quantum phase 
between the Neel and stripe phases. By contrast, for the 
s = 1 case, the gs energy curves of the Neel state for all 
values of J[ cross or meet the gs energy curves of the stripe 
state within the same physical region described above. 

Figure 2 shows our corresponding extrapolated results 
for the gs staggered magnetisation M. The quantum phase 
transition or critical point (Ma^.) marking the end of ei- 
ther the quantum Neel state or the quantum stripe state 
for a given value of J{ is first determined by calculating 
the order parameter M to obtain the value of J2 where 
M vanishes. However, as seen in fig. 2, there also oc- 
cur cases where the order parameters curves for the two 
states cross before their respective vanishing points. In 
such cases we take the crossing point to indicate the phase 
boundary between the quantum Neel and quantum stripe 
states. Thus, our definition of the quantum critical point 
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Fig. 2: Extrapolated CCM SUBn-^^ results for the gs staggered 
magnetisation, M, of the spin-1 Ji-J(-J2 model (with Ji — 1), 
for Ji = 0.2, 0.4, 0.6, 0, 7, 0.8, 0.9, 1.0. The SUBn-n results are 
extrapolated in the limit n — > 00 using the set n = {2, 4, 6, 8}. 



is the point where there is an apparent first-order phase 
transition between the two states or where the order pa- 
rameter vanishes, whichever occurs first. A fuller discus- 
sion of this former criterion and its relation to the stricter 
energy crossing criterion is given elsewhere [16]. 

We note particularly the result for this s = 1 model 
that the order parameter M curves for both the quan- 
tum Neel and stripe phases with the same value of J( 
go to zero smoothly at the same point, for all vahics of 
J[ < 0.66 ± 0.03. We emphasize that this cannot be acci- 
dental since it occurs for a large number of essentially in- 
dependent calculations over a wide parameter range. We 
also take this as further strong evidence for the accuracy 
and robustness of our extrapolation scheme. Thus, in this 
regime we have behaviour typical of a second-order phase 
transition between the quantum Neel and stripe phases. 
Furthermore, the transition occurs at a value of J2 very 
close to the classical transition point at J| = 0.5 J[. Con- 
versely, for values of J( > 0.66±0.03, the order parameters 
M of the two states meet at a finite value, as is typical of 
a first-order transition. 

Figure 3 shows the zero-temperature phase diagram of 
the spin-1 J\-J[-J2 model on the square lattice, as ob- 
tained from our extrapolated results for both the gs energy 
and the gs order parameter M. Unlike the spin- 1/2 case 
there is no sign at all of any intermediate disordered phase 
for any values of the parameters J{ or J2 (for Ji = 1). 
Hence, in this respect it is just like the classical case. How- 
ever, unlike the classical case, there now appears to be a 
quantum tricritical point (QTCP) at ,][ ss 0.66 ± 0.03, 
J2 « 0.35 ± 0.02, where a line of second-order phase tran- 
sitions between the quasi-classical Neel and stripe-ordered 
phases (for J{ < 0.66) meets a line of first-order phase 
transitions between the same two states (for J[ > 0.66). 
We note that the behaviour of both the order parameter 
(which goes to zero smoothly at the same point for both 
Neel and stripe phases below the QTCP, but which goes 
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Fig. 3: Extrapolated COM SUBn-n results for the gs phase 
diagram of the spin-1 Ji-J(-J2 model (with Ji = 1), showing 
a quantum tricritical point (QTCP). The SUBn^ results are 
extrapolated in the limit n — > oo using the set n = {2, 4, 6, 8}. 



to a nonzero value above it) and of the gs energy curves for 
the two phases (which meet smoothly with the same slope 
below the QTCP, but which cross with a discontinuity in 
slope above it) tell exactly the same story, as observed in 
other similar cases [16]. 

Thus, there is no evidence from our work for an interme- 
diate phase (for larger values of J[/Ji) for the s = 1 case, 
by contrast with the s = 1/2 case from our own previous 
results [17] (from which we found an intermediate phase 
without magnetic LRO for J[/Ji > 0.60 ±0.03) and those 
of other groups [24-26]. For the spin-1/2 case this interme- 
diate magnetically disordered phase was shown by us [17] 
to exist for the pure J1-J2 model (i.e., with J( = Ji) in 
the parameter range a J < a < for a = J2/J1, where 
al w 0.44 ± 0.01, « 0.59 ± 0.01, in full agreement with 
the accepted values. By contrast, for the s = 1 model pre- 
sented here we find instead a QTCP at J[/Ji ~ 0.66±0.03, 
J2/J1 « 0.35 ±0.02. 

For the case of the isotropic J1-J2 model lowest-order 
(or linear) spin-wave theory (LSWT) [1] predicts that 
quantum fluctuations can destabilize the classical GS with 
LRO even at large values of the spin quantum number 
s, for values of a = J2/J1 around 0.5. For the spin- 
1/2 case the range of values a J < a < a^, for which 
a magnetically-disordered phase occurs is predicted by 
LSWT to be given by « 0.38, ^ 0.52. For the corre- 
sponding spin-1 case LSWT predicts a much narrower, but 
still non-vanishing, strip of disordered intermediate phase 
with al ~ 0.47, f« 0.501. However, in an important pa- 
per, Igarashi [38] has shown explicitly for the spin-1 case, 
by going to higher-order terms in the 1/s power expansion 
of spin- wave theory (SWT), that no predictions based on 
LSWT (or SWT more generally) can be relied upon for 
values J2/J1 ^ 0.4 since the series seems to diverge in this 
region, with second-order terms becoming exceptionally 



large. Igarashi also showed that the higher-order correc- 
tion terms to LSWT act to make the Neel-ordcrcd phase 
more stable than LSWT would predict. We note too that 
Read and Sachdev [39], using a large- A/' expansion tech- 
nique based on symplectic Sp(A^) symmetry, which can it- 
self be regarded as akin to a 1/s expansion, have also found 
for the isotropic J1-J2 model an intermediate phase (with 
valence-bond solid order) for smaller values of s, which 
disappears for larger values of s where they find instead a 
first-order transition between the Neel and stripe phases. 
All of these results for the pure Ji- J2 model are in accord 
with ours. 

Naturally, one can also validly argue that what we have 
observed as a continuous (second-order) transition below 
the QTCP (i.e., for J(/ Ji < 0.66) might actually be a very 
weak first-order transition, which would thereby still com- 
ply with the Landau symmetry criterion of the standard 
Ginzburg-Landau theory of phase transitions and critical 
phenomena. Our completely independent sets of CCM cal- 
culations based on the two quasi-classical phases can never 
entirely exclude this possibility. However, our results from 
figs. 1 and 2 show clearly that the data below the QTCP 
are really only consistent with a transition which, if it is 
not second-order is at best very weakly first-order for all 
values of J^/Ji below the QTCP. In this context it is rele- 
vant to mention again that it has also been argued by oth- 
ers [12] that for the equivalent spin-1/2 model the phase 
transition between the Neel state and the intermediate 
paramagnetic state (which has been argued by those au- 
thors to be a valence-bond solid state) is also second-order 
and hence not described by standard Ginzburg-Landau 
critical theory. Again, for the spin-1/2 pure Ji J2 model 
the standard view is that the quantum phase transition 
between the striped and magnetically disordered interme- 
diate phases is first-order, and there is no discussion in the 
literature of deconfined quantum criticality for this tran- 
sition. One might argue, on similar grounds, that for our 
s = 1 model a first-order transition for the stripe phase 
might be more likely than for the Neel phase. We stress 
again, however, that our own results do indicate a direct 
second-order transition between these two phases below 
the QTCP. 

In a similar vein one might wonder too whether for the 

present spin-1 Ji J{ J2 model there might exist a nar- 
row strip of some intermediate phase, which could perhaps 
also act to reconcile our results with standard Ginzburg- 
Landau theory. Again, such a possibility cannot be ruled 
out with complete certainty by any numerical calculation 
such as ours. However, we have shown that our own ex- 
trapolation schemes are sufficiently robust and show suffi- 
cient internal consistency to rule out any but a very nar- 
row strip of an intermediate phase for < J{/Ji < 0.66. 
We estimate that the width of such a strip cannot ex- 
ceed by more than a factor of three or so that shown in 
fig. 3 from the data used in the present extrapolation. 
However, we note that in the limiting case of the spin-1 
ID chain (where J( — > 0, J2 ^ 0) the actual GS is the 
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Haldane gapped state [20] . Presumably this state should 
persist for small enough perturbations corresponding to 
small nonzero values of J[ and J2. The only other numer- 
ical study of the spin-1 Ji-J{-J2 model of which we arc 
aware [26,27] focused particular attention on this regime, 
and did indeed observe the continuation of the Haldane 
phase in a narrow strip in this regime. Our own results are 
not inconsistent with these observations, but our interest 
here lies more in the case of stronger interchain couplings 
where J{/Ji and J2/J1 are not confined to be small. How- 
ever, we note that other CCM calculations aimed specif- 
ically at this regime do, indeed, detect the Haldane gap. 
Thus, Zinke et al. [35] investigate the magnetic LRO of 
weakly coupled (quasi- ID) Heisenbcrg antifcrromagnetic 
chains for both the spin- 1/2 and spin-1 cases, using the 
CCM to calculate the staggered magnetisation and its de- 
pendence on the interchain coupling strength (J_l). They 
find that for the s = 1/2 case an infinitesimally small J± 
suffices to stabilize magnetic LRO, whereas for the s = 1 
case a nonzero (albeit small) J± is needed to establish 
LRO, in agreement with the results from other methods. 

Finally, in reaching our conclusions we have relied on 
two of the unique strengths of the CCM, namely its ability 
to deal with highly frustrated systems as easily as unfrus- 
trated ones, and its use from the outset of infinite lattices. 
There is no doubt that the many-body system studied here 
is highly non-trivial, and one cannot perhaps expect any 
single analysis or method to solve it completely. However, 
wc present for the first time new and interesting results 
using a method for which much previous work has shown 
its ability to describe quantum phase transitions reliably. 
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